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We study numerically the evolution of energy-level statistics as an integrability-breaking term 
is added to the XXZ Hamiltonian. For finite-length chains, physical properties exhibit a cross- 
over from behavior resulting from the Poisson level statistics characteristic of integrable models to 
behavior corresponding to the Wigner-Dyson statistics characteristic of the random-matrix theory 
used to describe chaotic systems. Different measures of the level statistics are observed to follow 
different crossover patterns. The range of numerically accessible system sizes is too small to establish 
with certainty the scaling with system size, but the evidence suggests that in a thermodynamically 
large system an infinitesimal integrability breaking would lead to Wigner-Dyson behavior. 

PACS numbers; 75.10.Pq, 75.40.Gb 



I. INTRODUCTION 

The conjecture that the statistical properties of energy 
levels of chaotic quantum systems may be described in 
terms of the theory of random matrices is widely accepted 
in various fields of physics i This however is not a uni- 
versal property of all complex interacting systems. One 
example to the contrary is provided by the class of the 
so-called integrable models,^ where the behavior of the 
system is completely described by a large (infinite in the 
thermodynamic limit) set of conserved quantities. One 
consequence is that the level-spacing distribution Pa{E) 
in the case of integrable models is the Poisson distribu- 
tion (A denotes mean level spacing). 



Pa{E) = -e 



-E/A 



(1) 



whereas in random-matrix theory the distribution takes 
the Wigner-Dyson form. 



PaIE) 



(2) 



where /3 = 1,2,4 correspond to orthogonal, unitary and 
symplectic ensembles respectively and^ bi — tt/2, ai — 
7r/4; &2 = 32/7r2 « 3.24, 03 = A/n; 64 = 262144/7297r3 « 
11.6, a4 = 64/97r « 2.26. 
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Other statistical properties (for example, the evolution 
of levels under an external perturbation"*) also differ for 
the two cases. One important class of external pertur- 
bations is the application of a voltage. The difference 
in response in this case leads to spectacular differences 
in transport properties of integrable and non-integrable 
models. Integrable models have been argued to have an 
infinite conductivity even at high temperatures, essen- 
tially because a typical level has a large response to a 
voltage, whereas non-integrable models have a finite con- 
ductivity because a typical level has a small response 

While these basic properties have been established for 
the two generic cases of integrable and non-integrable 
models, the crossover between these two limits as an 
integrability-breaking interaction is turned on has not to 
our knowledge been carefully studied, nor have the im- 
plications of the crossover for the finite-size conductivity 
of nearly integrable systems been determined. Two of 
us, with N. Andrei, presented a few numerical results 
in a paper mainly concerned with the charge transport 
of integrable systems.^ However, the significance and in- 
terpretation of these results was not clear. Song and 
Shepelyansky^ studied the effects of a random poten- 
tial on level statistics of 2D interacting Fermions and 
found evidence for a localization-delocalization transi- 
tion. However, in their case, the physics of the transi- 
tion is due to the disorder and thus is different from the 
situation in integrable models. Earlier work by DiStasio 
and Zotos^ noted a crossover between Poisson to Wigner- 
Dyson only for the low-energy part of the spectrum and 
did not address scaling with system size. 

In this paper we will fill these gaps by providing numer- 
ical results for finite-size chains with Hamiltonian given 
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by the (integrable) XXZ model plus an integrability- 
breaking perturbation SH. Our principal results are com- 
putations, for finite-length chains, of the crossover from 
behavior characteristic of Poisson to behavior character- 
istic of Wigner-Dyson statistics in various statistical mea- 
sures. These crossovers do not display an obvious univer- 
sality in the sense that different measures show different 
behavior depending on the XXZ asymmetry parameter 
and system size. 

Our computations are performed for finite-size sys- 
tems. An important issue is the behavior in the limit of 
thermodynamically large system. Extrapolation to the 
thermodynamic limit proves to be ambiguous for most 
of the measures we employ (namely, we cannot rule out 
a saturation of the crossover scales as functions of the 
system size for chains much longer than those considered 
in this study), but the data suggest that all the crossover 
scales vanish at infinite system size. 

The rest of the paper is organized as follows. First 
we discuss the model used in numerical calculations 
and in particular define numerically the value of the 
integrability-breaking parameter at which a gap appears 
in the spectrum. All further considerations will be de- 
voted to the gapless regime. Then we discuss the level- 
spacing distribution and the correlator of level veloci- 
ties. The latter is related to the parametric statistics 
of the system and also to its transport properties. For 
disordered systems, the correlator of level velocities was 
shown to correspond to the dimensionless conductance 
of the system, while if one restricts the analysis to pe- 
riodic boundary conditions only (see below) it coincides 
with the Drude weight. Discussion of the Drude weight 
concludes the paper. 

II. THE MODEL 

We study the effect of integrability breaking on the 
physical properties of a spin chain. The integrable model 
we consider is the XXZ chain defined on a A^-site ring 
with periodic boundary conditions in the presence of ex- 
ternal flux threading the ring: 



N 

+ ^ JiS'f S'f+i- (3) 

1=1 

(Alternatively, the flux can be gauged out to the bound- 
ary, resulting in twisted boundary conditions.) As is 
well-known, statistical properties of integrable models are 
governed by the Poisson distribution, Eq.lQJ. Transport 
properties of the model can also be inferred from studying 
the energy levels of the model, namely by their response 
to the flux <j). At zero temperature the behavior of the 



ground-state energy of the system under slow variations 
of the flux determines the Drude weight or the stiffnessiS 

Vs as 



Non-vanishing Vg signals ballistic transport in the sys- 
tem. For the XXZ model at T = this is the caseii for 
— 1 < Ji < 1, where excitations of the system are gapless. 
If |Ji| > 1, then the excitation spectrum of the model 
is gapped ji^ and T>s = 0. At finite temperatures the 
above expression for the stiffness can be generalized^*Sii 
to Vs = Di + D2, where 

_ A^ 1 d^Z 

vanishes in the thermodynamic limit, ^'^ and the rest is 
positive: 




In the gapless phase of the XXZ model it has been 
shown^' -'^ that ballistic transport persists to finite tem- 
peratures in the sense that D{N) — limr^oo D2{N)T > 
0. The infinite-temperature limit of this result implies 
that for a typical level dEn/dcj) ^ 1/Vn. At the antifer- 
romagnetic Heisenberg point Ji — I the model still has 
gapless excitations, but results of Fabricius and McCojii^ 
suggest that D2 vanishes (slowly) as the system size in- 
creases. Numerical results of Narozhny et al7 were con- 
sistent with this suggestion, but the limited range of sys- 
tem sizes attainable precluded a definite statement. 

Integrability breaking is introduced by adding the term 
with next-neighbor coupling 

N 

6H = Y,J2StSt^2 ■ (5) 

i=l 

This term should be contrasted to that considered by 
Eggerti^ insofar as it is explicitly not SU(2) invariant. 
However, away from the Heisenberg point the effect of 
the interaction Eq. |SJ) is similar to that of its SU(2)- 
invariant counterpart: (i) it breaks the integrability of 
the system, and (ii) for large enough values of J2 it causes 
the system to dimerize, so that the spectrum acquires a 
gap. The critical value of J2 at which the gap opens is 
of course different from the 0.24 found in Ref. [la- Our 
numerical estimatesi^ suggest a value Jj'^'* > 1.1 that is 
a weakly increasing function of both system size and Ji : 
for A = 18 and Ji — 0.2, for example, the gap appears 
to open at J2 = 1.13 ± 0.01. The gap opening limits the 
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range of values of J2 under consideration, as we are inter- 
ested only in properties of the gapless phase; indeed, the 
measures we consider presently begin showing different 
behavior for J2 > J2 • Similarly, the parameter region 
considered by Faas et alA belongs to the gapped regime, 
which accounts for certain differences in the behavior of 
the level statistics reported in Ref. and in the present 
paper. 

In this paper we study the eigenvalues of H ^ Hxxz + 
5H and their evolution under the change of (j) for the 
above models with varying Ji_2 and system size. For the 
non-integrable Hamiltonian Hxxz + 5H we use exact 
numerical diagonalization to construct the level-spacing 
distribution and level auto-correlation functions and to 
evaluate the stiffness I?2, Eq. Q). The use of exact nu- 
merical methods is motivated by the need to obtain the 
whole spectrum of the model in order (i) to study the sta- 
tistical properties of the spectrum and (ii) to study the 
stiffness Eq. (0J at infinite temperature. The drawback 
of the method is the limitation to small system sizes (we 
present results for chains of up to 20 sites). For finite 
system sizes we obtain a detailed characterization of the 
crossover. 
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FIG. 1: Typical crossover of the level-spacing distribution 
from Poisson (left solid curve) for J2 = to Wigner-Dyson 
(right solid curve) for a representative system. The plot is 
made for iV = 18, = 3, Ji = 0.2, and J2 = 0.1, 0.2, 0.5. For 
J2 = the numerical distribution agrees very closely with the 
exponential plotted. The Wigner-Dyson distribution shown 
is the theoretical curve for the orthogonal ensemble. 



III. LEVEL-SPACING DISTRIBUTION 

We begin with a brief discussion of the integrable case. 
The level-spacing distribution for J2 = is the Poisson 
distribution (shown in Fig^by the left solid curve). This 
illustrates the fact that the integrable system has so many 
conservation laws that levels essentially do not repel each 



other. To characterize transport properties of the system 
we show in Fig.Elthe quantity D{N) = D2{N)T at T ^ 
00 for different system sizes and different values of the 
integrable interaction Ji (dashed lines in Fig. [21) . I? is 
seen to be almost size-independent for the cases Ji < 
1, in agreement with previous work^i while a weak size 
dependence is evident in the Heisenberg case Ji = 1. 
Although this dependence appears to have a positive y- 
intercept, we believe that the system size in this study is 
still too small to make a definite statement regarding the 
behavior of the Heisenberg model in the thermodynamic 
limit 
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FIG. 2: The stiffness D{N) = limr^oo -D2(A^)r as a func- 
tion of inverse system size for different values of interaction 
parameters. Dashed lines show a naive extrapolation to the 
thermodynamic limit for the integrable system. Here we show 
the results for Ji — (squares), Ji = 0.2 (triangles) and 
Ji = 1 (circles). Empty symbols correspond to the smallest 
J2 = 0.01, which seems to affect strongly only the Ji = 
case. The behavior for J2 = 0.9 appears to be independent of 
Ji. 



We turn now to the case of broken integrability. As 
the integrability-breaking term Eq. (O is added to the 
ffamiltonian, energy levels immediately start to repel,— 
and as a consequence immediately ^'a(O) = so that 
the distribution acquires a peak. As illustrated in Fig.^ 
increasing J2 shifts the peak to the right until the distri- 
bution starts to look like the Wigner-Dyson distribution 
(shown in Fig. ^ by the right solid curve). At the same 
time the tail of the distribution changes from the expo- 
nential in Eq. to the (asymptotically) Gaussian tail 
of the Wigner-Dyson distribution. 

To quantify this crossover we show the evolution of 
the peak position and the characteristics of the tail with 
the change in J2 in Figs. 13 and 01 Both exhibit similar 
features, although the estimates for the crossover scales 
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extracted from the two are numerically different (see ta- 
ble U and insets in Figs. O and O. 




J. 

FIG. 3: Typical crossover of the peak position. The data 
correspond to A'^ = 20, = 4, Ji = 0.2, with the solid line 
a fit to the form atanh(J2/J2 ), J2 ~ 0.25. The dashed line 
indicates the peak position of the Wigner-Dyson distribution. 
The inset shows finite-size scaling of the crossover scale (the 
data points correspond to TV = 20, 5*^ = 4; iV = 20, = 3; 
Af = f8, S*" = 3; iV = 16, 5"" = 3; iV = f6, 5" = 2; iV = f4, 
= N = 14, = 2); the straight lines are guides to the 
eye, suggesting that J| vanishes for the infinite chain. 

As shown in Fig.|31 the peak of the distribution grows 
from zero to the value characteristic of the Wigner- 
Dyson distribution and then saturates. To estimate the 
crossover scale J2, we fit the data by the hyperbolic 
tangentiS of the form atanh(a;/xo) with xq approximat- 
ing J2. The inset shows the resulting values for as a 
function of the system size. As we noted beforeji^ we are 
restricting our attention to fixed values of the total spin 

. However, for the purposes of the finite-size scaling, 
it makes more sense to compare data with the fixed ratio 

/N. One way to see this is to recall that by means 
of the Jordan- Wigner transformation the spin chain can 
be mapped onto a system of spinless Fermions.^ In the 
Fermion language, 1/2 — /N corresponds to the filling 
fraction. Since it is not possible to keep the ratio /N 
exactly the same for all values of N used in this paper, 
we choose to present the data for two sectors of fixed 
that are closest to the chosen value of /N. Therefore 
the inset in Fig.|31shows two data points for the N other 
than N = 18 (we chose S'^/N = 1/6). The straight lines 
are just guides to the eye. 

To analyze the evolution of the tail, we approximate 
the intermediate distributions (see Fig. by 




Clearly, for Eq. Q a = 1 and 6 = 0, while for the orthog- 



onal ensemble, Eq. corresponds to a = and 6 = 7r/4. 
In Fig. 21 we show the evolution of b (the fact that plotted 
values never reach 7r/4 is an artifact of the calculation). 
Fitting the curve to a hyperbolic tangent, we can extract 
an estimate for the crossover value J2{N = 20) = 0.27. 
This value differs somewhat from the one extracted from 
the peak position (for the same values of N, S^, and Ji); 
see the inset in Fig.O The behavior of the tail character- 
istics with respect to changing system size exhibits the 
same trend as shown in the inset in Fig. |21 for the peak 
position: the characteristic scales tend to decay with in- 
creasing system size. The naive extrapolation of such a 
trend is consistent with a statement of vanishing as 
N — !■ Qo; however, the data are insufficient to prove it. 
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FIG. 4: Evolution of the parameter b from the tail of the level- 
spacing distribution. For the integrable case, 6 = 0. For large 
J2 it does not quite reach the Wigner-Dyson value 7r/4, but 
it still shows a clear crossover. The crossover scale estimated 
by fitting the data to a hyperbolic tangent is JI = 0.27. The 
plot is made for iV = 20, 5"= = 3, Ji = 0.2. 

Another way to quantify the evolution of the level- 
spacing distribution shown in Fig. ^ is to consider 
cumulants^ Their unbiased estimators (the Fisher 
statisticsSi kn) are easily computed. For our normalized 
level spacings, the first cumulant (which is equal to the 
mean) is unity. 

The cumulants of a distribution characterize its width 
(second cumulant, or variance) and shape;^^ beyond per- 
haps the fifth, numerical cumulants become too sensitive 
to outliers to be of much use. A study of the cumulants 
of a distribution is qualitatively similar to our foregoing 
study of the tails, but it turns out to be simpler numeri- 
cally. In Fig. [S] we show the unbiased variance estimate, 
A:2, as a function of J2 for system size N = 20, Ji = 0.2, 
S"^ = 3, momentum k = 0. The theoretical limits should 
be 1 for the Poisson distribution and 7r/4 — 1 = 0.2732 
(the bottom of the scale) for Wigner-Dyson. The fact 
that the data points deviate from these ideal values is 
an artifact of finite sampling. Fitting the curve to a hy- 
perbolic tangent, we estimate the crossover scale, shown 
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in the inset as a function of size (for the same sequence 
of quantum numbers as in Fig.|2J). The inset also shows 
the crossover scale for the fourth cumulant. As before, 
the scaling suggests (but does not establish) that the 
crossover scale J2 associated with either cumulant should 
vanish in the limit of infinite size. It is not possible to 
determine whether the different measures, k2 and (we 
also looked at k^), scale in the same way or differently 
with system size. 
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FIG, 5: Variance estimate (x), fe, as a function of J2 for 
system size A'^ — 20. (Ji = 0.2, = 3.) The crossover scale 
is estimated as the turning point (o) in a tanh fit (curve). The 
inset illustrates the finite-size scaling for cumulants k2 (solid 
lines and x) and ki (dotted lines and <>) for the same set of 
sizes and sectors as in Figs. |3 and |7| appropriate for 1/3 
filling. The JI associated with either cumulant may vanish 
in the limit of infinite size. 

The second cumulant k2 plays a role similar to that of 
the parameter 77 used in Ref. \§ to estimate the overall 
"proximity" of the observed distribution to either the 
Poisson or the Wigner-Dyson limit. In that sense, Fig.jS] 
shows behavior similar to that found in Ref. \§, although 
the physics of the evolution of levels is quite different in 
our case. 



IV. ELEMENTS OF THE PARAMETRIC 
STATISTICS 

More information about the crossover to the chaotic 
behavior described by the Wigner-Dyson statistics can be 
extracted from the study of the autocorrelation functions. 
Here we will discuss the autocorrelation of level velocities 

where the angular brackets indicate averaging over a set 
of levels and fluxes 1^ 
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FIG. 6: Autocorrelation of level velocities for various values 
of integrability breaking. Only a short range of flux is shown; 
in all cases, the sum rule (the integral of C{(j}) over all fluxes 
vanishing) is satisfied numerically. For J2 > 0, there is a point 
of inflection at non-zero flux. However, the integrable model 
has no initial inflection. 



A typical form of the autocorrelation function Eq. © 
is shown in Fig. IHl For large values of J2, this form re- 
sembles the universal correlator characteristic of chaotic 
systems J. However, in the crossover region, C{<j>) devi- 
ates from the universal form in a rather complex fashion, 
which makes a quantitative analysis of the crossover dif- 
ficult. Therefore we focus on two particular features of 
the curve, the turning point and C(0). We note a fea- 
ture of the autocorrelation curve shown in Fig. that 
appears only as the integrability breaking is introduced: 
all curves for J2 > have a non-zero point of inflec- 
tion as the autocorrelation decreases with the increase 
of J2 , but the autocorrelation function of the integrable 
system does not have such an inflection point. Conse- 
quently, the behavior near zero flux changes from linear 
(C(0) — C(0) cx —(p) for the integrable case to quadratic 
for J2 > 0. 

The autocorrelation function of level velocities at zero 
flux difference C(0) is somewhat similar to the stiffness 
Eq. 10} , the differences being that C(0) is also averaged 
over a set of fluxes, does not contain the extra factors of 
temperature and system size, and corresponds to a single 
sector of fixed S^. However, in chaotic systems it is C(0) 
that can be related to transport There it was argued to 
correspond to the dimensionless conductance. 

In Fig. [3 we show the behavior of C(0) as a function 
of J2- Clearly, for finite systems, C(0) exhibits a well- 
defined crossover. For the dataset presented in Fig. [3 
{N = 20, Ji = 0.2, 5'^ = 3), the correlator C(0) decays 
as a w 5/2-power law after J2 exceeds the value Jj*^" = 
0.097± 0.004 (defined in Fig.[7|as a crossing point of the 
above power law — the straight line in the log-log scale 
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reduction in the stiffness for finite chains. The extrapola- 
tion to infinite size is also uncertain. It should be noted, 
however, that if one compares the behavior of two inte- 
grable cases, Ji — and Ji > 0, then a similar picture 
arises (compare for example the two top dashed lines in 
Fig. [2). This has to do with the fact that the spin chain at 
Ji = can be mapped (by means of the Jordan- Wigner 
transformation) onto a system of free spinless Fermions 
and as such possesses more symmetries than even the 
integrable (but interacting) XXZ model. 

(iv) When the integrability-breaking parameter is not 
small, the stiffness decays sharply (in fact, if we were to 
show a log-log plot, faster than any power law) with sys- 
tem size and clearly extrapolates to zero in the thermody- 
namic limit. This behavior is qualitatively independent 
of Ji. 



FIG. 7: C(0) as a function of J2 for TV = 20, Ji = 0.2, S"" = 3. 
The dotted line corresponds to the value at J2 = 0. The inset 
shows the finite-size scaling of the crossovGr sc&ls, J2 ^ i 
approximately fixed ratio S^N = 1/6 (for iV = 14, 16, or 20, 
we pick the two closest values of 5*^ for these sizes). 

— with the value at J2 = 0). 

The inset shows the crossover scale as a function of 
the size (in the same manner as the crossover scale ex- 
tracted from the peak of the level-spacing distribution). 
The behavior is very similar to that in the inset in Fig. O 
(although numerical values of the crossover scales differ 
in the two cases). Both would be consistent with the 
statement that J| ^ as ^ 00; however such a con- 
clusion cannot be ascertained on the basis of the data 
available. 



V. SPIN STIFFNESS 

Now we discuss the effect of the integrability breaking 
Eq. (O on the stiffness D2 ■ In Fig. |21 we show three 
sets of data corresponding to three different values of the 
XXZ anisotropy parameter Ji. The data illustrate the 
following tendencies: 

(i) For < Ji < 1 (represented by Ji = 0.2; similar be- 
havior is observed for other values) the data clearly show 
that very small integrability breaking (characterized by 
J2 = O.OI) has little effect on the stiffness of the finite 
chains (which is to be expected). Moreover, the extrapo- 
lation to infinite size seems to result in a finite value for 
the stiffness in a manner similar to that of the integrable 
model. 

(ii) For the Heisenberg model, small integrability 
breaking again does not have a pronounced effect; how- 
ever, in this case (even though the extrapolation indicates 
a finite value for the thermodynamic limit), one cannot 
be certain of the behavior of the infinite chain. 

(iii) For Ji = the situation is different: even a very 
small amount of integrability breaking leads to a sharp 
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FIG. 8: D{N) as a function of J2 (for Af = 18). The main 
plot and the left inset show the same data in linear and log- 
log scale (meant to illustrate the crossover to a power-law 
decay of the stiffness at large J2 for smaller values of Ji) for 
Ji = (squares), Ji = 0.2 (circles), and Ji = 1 (diamonds). 
The horizontal lines in the inset indicate the values of D{N) 
at J2 = for Ji — 0.2 (the upper line) and Ji = 1. The 
second inset shows the size dependence of the crossover scale 
for Ji = 0, 0.2, and 0.4 (triangles). The Heisenberg case 
Ji = 1 does not show a clear crossover. 

Thus the data suggest that for finite chains there exists 
a "critical" value of J2 (smaller than the point of the gap 
opening) beyond which the stiffness tends to vanish. This 
value is not universal, in the sense that it depends on Ji. 
This situation is illustrated in Fig. |S1 where we show the 
dependence of the stiffness on J2 for three values of Ji 
and the fixed system size N = 18. For large values of 
J2, all three curves saturate to zero (although the one 
with Ji ~ 1 does so faster). For Ji = 0.2 the effect of 
small J2 is rather weak, and the curve exhibits a clear 
crossover. For Ji — the crossover also appears; how- 
ever, the value 0.08 to which the curve tends as J2 ^ 
is much smaller than the value at J2 = exactly (which 
is 0.25 and is thus outside the frame of the plot). The 
crossover is illustrated in the left inset. The right in- 
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set shows the change of the crossover scale with system 
size for Ji — 0, 0.2, and 0.4. The behavior at Ji = 1 is 
quite different and in particular does not show an obvious 
crossover (and is thus not represented in the right inset). 
This behavior might be related to the conjecturc^^ that 
at the Heisenberg point the stiffness vanishes in the ther- 
modynamic limit even without the integrability breaking. 
Alternatively, this can reflect the fact that the Heisenberg 
model is characterized by logarithmic correlations^ and 
thus the small systems considered in this paper are not 
representative. 

VI. DISCUSSION 

Prior to performing the calculations one could have 
had two conflicting expectations for the behavior of the 
non-intcgrable system: (i) as soon as the integrability is 
broken the system becomes chaotic and as a result shows 
diffusive transport; (ii) there exists a "critical" mag- 
nitude of the integrability breaking that separates the 
chaotic regime from the one that retains some features of 
the integrable model, in particular, ballistic transport. 

The latter picture has an analogy in the localization 
problem in disordered conductorsi^^ The states of an in- 
tegrable model can be visualized as well-defined local- 
ized points in the multidimensional space of the integrals 
of motion characteristic of the model. These points are 
well separated due to the quantization of the values of 
the integrals of motion. Consider now the effect of an 
infinitesimally small integrability breaking. One can cer- 
tainly expect the points to spread out into fuzzy spots, 
but at the same time one might argue that unless the in- 
tegrability breaking is strong enough, these spots do not 
overlap. In this regime the system retains some memory 
of the fact that it was indeed integrable before the extra 
interaction was turned on. When the integrability break- 
ing is so strong that the spots overlap into a continuum, 
the system becomes fully chaotic. 

The numerical analysis presented in this paper seems 
roughly consistent with the second possibility for fi- 
nite chains: a small integrability-breaking term leads 
to behavior that is close to that of the integrable sys- 
tem. Quantities related to transport, the stiffness D2 
and the "conductance", C(0), exhibit a reasonably rapid 
crossover as functions of the strength of the integrability- 
breaking interaction. The crossover behaviors seem to 
be different for different quantities. Table for exam- 
ple, illustrates some of this variability for the example 
of Ji = 0.2, intermediate between the non-interacting 
model and the Heisenberg point. Similarly, the quan- 
titative characteristics of the level-spacing distribution 
(namely the peak position and the tail parameter, see 
Figs. 121 and 21 or cumulants, see Fig. O exhibit similar 



crossovers (the corresponding scales are also included in 
table P). 

The one exception to this picture is the number of 
degenerate levels in the systemiSi represented by Pa(0). 

TABLE I: The system crosses over from integrable to fully 
chaotic behavior with different crossover scales, depending 
on what is being measured. Furthermore, the crossover scales 
themselves scale differently with system size. We calculate 
crossover scales associated with peak position, tail crossover 
(from exponential to Gaussian), mean squared level veloc- 
ity Co, the fourth cumulant k4 of the level-spacing distribu- 
tion, and conductance D2T. In this example, Ji = 0.2. The 
crossover Jj^^"^ is calculated for the entire spectrum, while 
all the others are calculated for the S'' — 3 sector. Entries of 
— could not be extracted from the data because of numerical 
uncertainty. Rough error estimates for the least significant 
digit are provided where available. 



N 


T*pcak 


T*tail 


T*Co 


"^2 


r*-D2T 


20 


0.19 


0.27 


0.097(4) 


0.091 


0.079(3) 


18 


0.25 


0.43 


0.15(1) 


0.17 


0.107(3) 


16 


0.34 


0.49 


0.20(2) 


0.24 


0.145(2) 


14 


0.38 




0.28(2) 




0.178(5) 



This measure exhibits a jump as infinitesimally small (nu- 
merically meaning of the order of the computer precision) 
J2 is introduced (namely -Pa(O) vanishes, as illustrated 
in Fig.jTJ. 

Conclusions for the thermodynamic limit are harder to 
draw from our data. The variation of crossover scale with 
system size indicated in the insets to the different figures 
suggests that the crossover scale vanishes in the limit of 
infinite size system rather than saturating at non-zero 
values for , but the limited range of sizes available to 
us, along with the absence of a theoretically justified ex- 
trapolation to the thermodynamic limit, precludes a def- 
inite statement. Constructing a theory of the approach 
to the infinite-size limit of chains with weak integrability 
breaking remains an important open problem. 
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